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J.  A.  NOVOTNY,  E.  C.  PARKER,  S.  S.  SURVANSHI,  G.  W.  ALBIN,  AND  L.  D.  HOMER 
Naval  Medical  Research  Institute,  Bethesda,  Maryland  20889-5055 


Novotny,  J.  A.,  E.  C.  Parker,  S.  S.  Survanshi,  G.  W.  Al- 
BIN,  and  L.  D.  Homer.  Contribution  of  tissue  lipid  to  long  xenon 
residence  times  in  muscle.  J.  Appl.  Physiol.  74(5):  2127-2134, 
1993. — Experiments  demonstrate  that  the  mean  residence 
time  of  an  inert  gas  in  tissue  is  longer  than  that  predicted  by  a 
single-compartment  model  of  gas  exchange.  Also  the  relative 
dispersion  (RD,  the  standard  deviation  of  residence  times  di¬ 
vided  by  the  mean)  is  1  according  to  this  model,  but  RDs  in  real 
tissues  are  closer  to  2,  suggesting  that  a  multiple-compartment 
model  might  be  more  accurate.  The  residence  time  of  a  gas  is 
proportional  to  its  solubility  in  the  tissue.  Although  the  noble 
gases  in  particular  are  10  times  more  soluble  in  lipid  than  in 
nonlipid  tissues,  models  of  gas  exchange  generally  do  not  incor¬ 
porate  measurements  of  the  lipid  in  tissue,  which  may  lead  to 
error  in  the  predicted  gas  residence  times.  Could  a  multiple- 
compartment  model  that  accounts  for  the  lipid  in  tissue  more 
accurately  predict  the  mean  and  RD  of  gas  residence  times?  In 
this  study,  we  determined  the  mean  and  RD  of  Xe  residence 
times  in  intact  and  surgically  isolated  muscles  in  a  canine 
model.  We  then  determined  the  lipid  content  and  the  perfusion 
heterogeneity  in  each  tissue,  and  we  used  these  measurements 
with  a  multiple-compartment  model  of  gas  exchange  to  predict 
the  longest  physiologically  plausible  Xe  residence  times.  Even 
so,  we  found  the  observed  Xe  mean  residence  times  to  be  twice 
as  long  as  those  predicted  by  the  model.  However,  the  predicted 
RDs  were  considerably  larger  than  the  observed  RDs.  We  con¬ 
clude  that  lipid  alone  cannot  account  for  the  residence  times  of 
Xe  in  tissue  and  that  a  multiple-compartment  model  is  not  an 
accurate  representation  of  inert  gas  exchange  in  tissue.  A  more 
accurate  model  will  need  to  account  for  other  features  in  the 
tissue,  such  as  diffusion  of  gas  between  vessels  and  structures 
with  heterogeneous  solubility. 

inert  gas;  mathematical  model;  decompression  modeling 


the  simplest  model  of  the  exchange  of  dissolved  inert 
gases  between  blood  and  tissue  considers  the  tissue  as  a 
single  well-mixed  compartment.  However,  empirical 
studies  on  inert  gas  exchange  find  that  the  mean  resi¬ 
dence  time  observed  in  the  tissue  is  nearly  always  longer 
than  that  predicted  by  the  model  (2,  7,  9).  According  to 
the  single-compartment  model,  the  gas  uptake  and  elimi¬ 
nation  curve  is  a  single-exponential  function,  and  the 
mean  residence  time  of  the  gas  is  equal  to  X/F,  where  the 
tissue-blood  partition  coefficient  (X)  is  the  ratio  of  the 
average  solubility  of  the  gas  in  the  tissue  to  its  solubility 
in  the  blood  and  F  is  the  blood  flow  rate  per  unit  volume. 
The  relative  dispersion  (RD,  the  standard  deviation  of 
gas  residence  times  divided  by  the  mean)  of  a  single-ex¬ 
ponential  process  is  1,  but  RDs  observed  in  tissues  are 
closer  to  2  ( 12, 15),  indicating  a  wider  distribution  of  resi¬ 


dence  times.  A  multiple-compartment  model  can  ac¬ 
count  for  RDs  >1  and  so  may  be  a  more  accurate  repre¬ 
sentation  of  inert  gas  exchange  in  tissue.  However,  we 
previously  found  that  such  a  model  can  account  for  only 
a  small  fraction  of  the  observed  RD  when  only  the  hetero¬ 
geneous  blood  flow  in  the  tissue  and  an  assumed  con¬ 
stant  value  for  X  are  considered  (6).  Is  another  model 
required,  or  is  the  multiple-compartment  model  in  error 
because  of  incomplete  information  about  the  parameters 
F  and  X? 

One  previous  study  on  Xe  exchange  in  skeletal  muscle 
(15)  attributed  the  unexpectedly  long  Xe  elimination 
times  to  bits  of  lipid  that  were  observed  on  visual  inspec¬ 
tion  of  the  tissue.  The  residence  time  of  an  inert  gas  is 
proportional  to  its  solubility  in  the  tissue,  and  the  noble 
gases  in  particular  are  ~  10-fold  more  soluble  in  lipid 
than  in  nonlipid  tissues  (17).  The  lipid  content  of  muscle 
tissue  is  variable  (4,  8),  so  tissues  with  above-average 
amounts  of  lipid  would  have  longer  gas  residence  times 
than  would  be  predicted  considering  only  the  average 
published  values  of  gas  solubility  in  tissue.  How  much 
lipid  would  be  necessary  to  account  for  the  Xe  residence 
times  observed  in  tissues? 

In  this  study,  we  tested  the  hypothesis  that  the  lipid 
present  in  a  sample  of  normal  canine  muscle  and  subcuta¬ 
neous  tissues  can  account  for  the  mean  residence  time  of 
Xe  gas.  We  substituted  values  for  X  into  a  multiple-com¬ 
partment  model  on  the  basis  of  the  measured  lipid  con¬ 
tent  of  the  tissue  and  biased  the  model  to  estimate  the 
longest  plausible  mean  residence  times.  We  also  esti¬ 
mated  the  effect  of  the  lipid  on  the  RD  of  residence  times 
in  the  presence  of  heterogeneous  tissue  blood  flow  rates. 

METHODS 

We  studied  eight  dogs  ( Canis  familiar  is,  1-2  yr  old, 
male,  13-26  kg)  anesthetized  with  pentobarbital  sodium 
and  mechanically  ventilated.  The  methods  were  similar 
to  those  described  previously  (13).  Only  modifications  to 
those  techniques  are  described  here. 

The  experimental  protocol  for  this  project  was  re¬ 
viewed  by  the  institutional  Animal  Care  and  Use  Com¬ 
mittee  and  certified  as  conforming  to  guidelines  in  Insti¬ 
tute  for  Laboratory  Resources,  National  Research  Coun¬ 
cil  “Guide  for  the  Care  and  Use  of  Laboratory  Animals” 
(DHHS  publ.  86-33). 

Xe  Administration 

The  dogs  breathed  from  a  closed-circuit  ventilator.  At 
the  start  of  the  experiment,  ,33Xe  gas  was  injected  into 

°1 


2128 


LIPID  AND  Xe  EXCHANGE  IN  MUSCLE 


the  ventilator  to  establish  a  step  increase  in  the  concen¬ 
tration  of  l33Xe  in  the  inspired  air,  which  was  then  main¬ 
tained  constant  for  2.5  h  of  Xe  uptake  by  continuously 
adding  t33Xe  and  monitoring  the  radioactivity  of  the  in¬ 
spired  limb  of  the  ventilator  circuit.  At  the  end  of  the  Xe 
uptake  period,  the  ventilator  was  suddenly  switched  to 
open-circuit  breathing  for  2.5  h  of  Xe  elimination. 

Tissue  Site  Preparation 

The  skin  and  subcutaneous  tissues  may  contain  consid¬ 
erably  more  lipid  than  the  underlying  muscle  tissue.  We 
studied  both  intact  hindleg  muscles  and  muscles  stripped 
of  the  overlying  nonmuscle  tissues  to  increase  the  range 
of  lipid  contents  of  our  sampled  tissues.  We  randomly 
selected  either  the  right  or  the  left  calf  muscle  for  surgi¬ 
cal  isolation.  The  skin,  subcutaneous  fat,  and  connective 
tissues  were  lifted  by  blunt  dissection  and  retracted,  and 
the  isolated  muscle  was  covered  in  0.6-mm-thick  plastic 
wrap  (Saran)  and  suffused  with  0.9%  saline  at  37°C.  If  we 
assume  that  the  plastic  wrap  (product  code  06042,  Dow 
Plastics,  Dow  Chemical)  is  no  more  permeable  to  Xe 
than  to  N2,  we  can  expect  <1 /800th  of  the  Xe  gas  deliv¬ 
ered  to  the  tissue  in  blood  to  be  lost  through  this  surface, 
a  small  loss  that  would  not  be  reliably  detected  with  our 
system.  The  other  leg  was  left  intact.  The  surface  temper¬ 
ature  of  the  isolated  muscle  and  the  subcutaneous  tem¬ 
perature  of  the  intact  side  (via  a  3-mm  incision  in  the 
skin)  were  continuously  monitored  with  2- mm  probes 
(Digitec,  Dayton,  OH).  Similar  temperatures  were  main¬ 
tained  on  both  sides. 

Xe  Detection 

Using  9-mm-diam  cylindrical  detectors  containing 
CdTe  semiconductors  (RMD,  Watertown,  MA),  we  re¬ 
corded  the  radioactivity  of  133Xe  as  a  function  of  time. 
Using  8-mm-thick  barrel-shaped  lead  collimators,  we 
shielded  the  detectors.  We  aligned  the  collimators  vi¬ 
sually  over  the  intact  and  the  isolated  muscle  prepara¬ 
tions,  excluding  the  bone  from  the  field  of  view  of  the 
detector.  We  simultaneously  recorded  the  background  ra¬ 
dioactivity  with  a  detector  positioned  midway  between 
the  two  detectors  over  the  tissues.  The  background  radio¬ 
activity  is  time  dependent,  reaching  ~500  cpm  at  the 
peak  of  an  experiment  compared  with  a  typical  tissue 
signal  of  10,000  cpm.  To  simulate  the  secondary  radia¬ 
tion  generated  by  the  ambient  radioactivity  and  the  de¬ 
tected  tissues  (15),  beneath  the  background  detector  we 
positioned  a  water-filled  plastic  cone  with  the  same  di¬ 
mensions  as  the  detector  field  of  view. 

We  subtracted  the  background  at  each  time  point  from 
the  tissue  data.  To  smooth  the  fluctuations  in  difference 
of  two  random  variables,  we  first  fitted  curves  to  the 
background  data  and  took  the  net  tissue  signal  as  the 
difference  between  the  raw  tissue  data  and  the  smoothed 
background. 

Measuring  the  Tissue  Layer  Thickness 

We  measured  the  thicknesses  of  the  muscle  and  non¬ 
muscle  tissues  of  the  intact  side  at  the  end  of  each  experi¬ 
ment.  This  information  is  important,  because  the  skin 
and  subcutaneous  tissues  are  closer  to  the  detector  and 


therefore  contribute  more  efficiently  to  the  detected  sig¬ 
nal  than  the  more  distant  muscle  tissue.  The  digital  cali¬ 
pers  had  a  resolution  of  0.01  mm  (Ultra-cal  II,  Fowler, 
Westmont,  IL),  which  resulted  in  a  standard  deviation  of 
<1  mm  on  repeated  measurements.  The  appendix  ex¬ 
plains  how  we  used  these  measurements  to  estimate  the 
contribution  of  the  nonmuscle  tissue  to  the  total  detected 
Xe  radioactivity  of  the  intact  side. 

Measuring  the  Tissue  Blood  Flow 

To  determine  an  average  blood  flow  for  each  tissue 
piece  during  the  5-h  Xe  exchange  period,  we  used  radioac¬ 
tive  microspheres  to  make  four  equally  spaced  measure¬ 
ments.  The  first  was  at  the  beginning  of  Xe  uptake,  with 
two  more  at  100-min  intervals,  and  the  last  was  at  the 
end  of  the  Xe  elimination  period.  We  used  15-pm  micro¬ 
spheres  labeled  with  51Cr,  95Nb,  103Ru,  or  141Ce  (NEN, 
Dupont,  Wilmington,  DE)  and  injected  them  in  aliquots 
of  2  X  105  microspheres/kg  in  10  ml  of  0.9%  saline  over  15 
s  through  a  catheter  placed  in  the  left  ventricle  via  the 
left  common  carotid  artery.  Reference  blood  samples 
were  simultaneously  withdrawn  from  a  brachial  artery 
catheter  at  2  ml/min  for  2  min  and  accompanied  by  an 
equal-volume  infusion  of  Dextran-70  into  a  pulmonary 
artery  catheter. 

At  the  end  of  the  experiment,  using  a  scalpel,  we  re¬ 
moved  the  tissues  in  the  field  of  view  of  the  Xe  detectors. 
We  then  manually  separated  the  tissues  of  the  intact  side 
into  muscle  and  nonmuscle  portions.  The  average 
weights  of  the  muscle  and  the  nonmuscle  tissues  were  30 
and  10  g,  respectively.  We  then  cut  the  tissues  into  1-g 
pieces  with  a  scalpel,  weighed  them,  and  placed  them  in 
plastic  tubes  for  counting  in  a  four-channel  scintillation 
detector  (LKB  Wallac,  Gaithersburg,  MD),  which  en¬ 
abled  the  estimation  of  the  blood  flow  rates  at  the  four 
times. 

The  average  tissue  blood  flow  over  the  four  measure¬ 
ments  for  the  ;th  piece  (F;,  ml  •  min  1  •  ml'1)  was  calcu¬ 
lated  as 

F,  =  l/4  7^-W  ( 1 ) 

k=l 

where  R;*  is  the  radioactivity  for  the  fcth  measurement  of 
the  yth  tissue  piece,  Rp*  is  the  radioactivity  of  the  arterial 
reference  sample,  and  W  is  the  reference  withdrawal  rate 
(ml/min).  The  volume  of  the  tissue  was  calculated  as 
described  below. 

Measuring  the  Tissue  Lipid 

After  counting  the  microsphere  radioactivity,  we 
pooled  the  tissue  pieces  into  their  three  original  group¬ 
ings  (muscle  from  the  isolated  side,  muscle  from  the  in¬ 
tact  side,  and  nonmuscle  from  the  intact  side)  and  ex¬ 
tracted  the  lipid.  The  limitations  of  the  lipid  extraction 
method  (1)  prevented  measurement  of  the  lipid  content 
of  the  individual  1-g  pieces  of  tissue.  The  lipid  contents  of 
our  tissue  samples  ranged  from  0.5  to  7.5  g.  We  tested  the 
accuracy  of  the  method  by  extracting  amounts  of  corn  oil 
in  the  same  range  as  our  tissue  lipid  values.  We  found 
that  the  oil  was  accurately  recovered  in  this  range,  with  a 
coefficient  of  variation  of  4%  on  replicate  samples.  Using 
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densities  of  0.9  for  the  lipid  and  1.1  for  the  residual  non¬ 
lipid  material  (11),  we  estimated  the  volumes  of  the  tis¬ 
sue  samples  according  to  the  lipid  fractions  and  the  mea¬ 
sured  weights. 


Analysis 

We  compared  the  observed  means  and  RDs  of  Xe  resi¬ 
dence  times  with  the  predictions  of  a  model  that  accounts 
for  the  lipid  and  the  blood  flow  of  the  tissue. 

Observed  Xe  kinetics.  The  Xe  radioactivity-time  curves 
less  the  smoothed  background  were  fitted  with  sums  of 
exponentials,  as  described  previously  (12).  The  mean  and 
the  RD  of  residence  times  were  then  calculated  from  the 
parameters  of  the  least-squares  fit. 

Predicted  Xe  kinetics.  The  model  includes  variables  to 
represent  both  the  lipid  content  and  the  blood  flow  heter¬ 
ogeneity  of  the  tissue.  For  the  intact  side,  we  estimated 
the  mean  and  the  RD  of  Xe  residence  times  as  averages 
of  the  muscle  and  nonmuscle  tissues  weighted  according 
to  their  volumes  and  relative  positions  beneath  the  de¬ 
tector  (appendix).  The  details  of  the  model  are  outlined 
in  Eqs.  2-13. 

For  each  experiment,  the  tissue-sampling  method 
yields  several  1-g  pieces  of  muscle  and,  on  the  intact  side, 
muscle  and  nonmuscle  tissue.  Each  1-g  piece  was  consid¬ 
ered  as  a  well-mixed  compartment  for  the  exchange  of 
Xe  between  blood  and  tissue.  The  composite  of  tissue  in 
the  detector  field  of  view  was  then  modeled  as  a  collec¬ 
tion  of  n  such  compartments.  The  yth  piece  has  an  expo¬ 
nential  time  constant  rj  (minutes)  that  depends  on  the 
partition  coefficient  and  the  blood  flow  per  unit  volume, 

F ) 


The  predicted  Xe  transfer  function  for  the  ith  tissue  is 


h(£),  =  Z  l^-jo/r,-)  exp  i-t/tj)  (3) 

with  i=  1-3  representing  the  muscle  tissue  of  the  isolated 
side  and  the  muscle  and  nonmuscle  tissues  of  the  intact 
side,  respectively.  FT  is  the  total  of  the  F;. 

The  predicted  Xe  mean  residence  time,  the  first  mo¬ 
ment  about  the  origin  (EM;)>  is  the  flow-weighted  aver¬ 
age  of  the  tissue  pieces 


EM, 


,-i  l  Ft 


( WF,>  + 


\Ft / 


(AnlV/F,)  (4) 


where  XNL  and  XL  are  the  tissue-blood  partition  coeffi¬ 
cients  for  the  nonlipid  and  the  lipid  compartments,  re¬ 
spectively,  and  nL  is  the  number  of  the  tissue  pieces  con¬ 
sidered  as  lipH  Although  our  lipid  extraction  method  did 
not  determine  the  lipid  contents  of  the  individual  pieces, 
the  mean  residence  time  depends  only  on  the  total  lipid 
content  and  the  total  blood  flow,  as  explained  below. 

The  tissues  we  harvested  from  the  animals  included 
muscle  tissue  and  skin  and  subcutaneous  nonmuscle  tis¬ 
sue,  for  which  the  average  published  tissue-blood  parti¬ 
tion  coefficients  for  Xe  are  0.68  (17)  and  0.87  (10),  respec¬ 
tively.  Our  sampled  tissues  had  a  wide  range  of  lipid  con¬ 
tents.  To  account  for  this  variability,  we  estimate  A  for 


the  tissues  in  each  experiment  according  to  their  mea¬ 
sured  lipid  contents.  We  assigned  these  published  values 
of  X  to  the  tissues  having  the  lowest  lipid  contents  of 
those  sampled  in  all  the  experiments.  We  then  accounted 
for  the  higher  lipid  contents  of  the  tissues  in  all  the  other 
experiments  by  assigning  XL  =  9.8  (17)  to  the  nL  pieces 
corresponding  to  the  volume  fraction  of  lipid  in  Eq.  4. 
This  ensures  that  all  the  sampled  tissues  have  average  X 
values  at  least  as  large  as  the  published  values  and  that 
the  model  estimates  mean  residence  times  as  long  as 
physiologically  plausible. 

How  does  the  predicted  mean  residence  time  depend 
on  the  spatial  arrangement  of  the  lipid  and  the  blood  flow 
in  the  tissue?  Summing  the  volumes  to  total  volume  ( VT) 
and  rewriting  Eq.  4  in  terms  of  the  tissue  lipid  fraction 
(w)  show  that  the  mean  residence  time  depends  only  on 
the  ratio  of  VT  to  Fx  and  the  weighted  average  of  the  XL 
and  XNL  but  not  on  their  location  in  the  tissue 

EM.  =  I1  [wXL  +  (1  -  U,)XNL]  (5) 

T  T 

We  can  use  this  expression  to  estimate  the  tissue-blood 
partition  coefficient  for  each  tissue  as  the  weighted  aver¬ 
age  of  the  lipid  and  nonlipid  portions 

X.  =  u)XL  +  (1  —  w)XNL  (6) 

where  X.  is  the  average  X. 

The  predicted  second  moment  about  the  origin  (EM]) 
of  the  distribution  of  Xe  residence  times  is 


EML  =  2I 
/■“  1 


(AlV,./F/  +  2 
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which  is  used  to  calculate  the  variance  and  the  RD. 

The  second  moment,  and  so  the  RD,  depends  on  the 
distribution  of  lipid  and  blood  flow  in  the  tissue.  This  can 
be  shown  by  simplifying  Eq.  7  to 


EML  =  2^Z  V,.T;  +  2-^  s 

rT;«l  rT  j  =  nt+l 


V, 


(8) 


Because  the  tissue  piece  volumes  (Vy)  are  similar  in  size, 
we  can  replace  them  with  a  standard  volume  (V9)  to  sim¬ 
plify  Eq.  8  to 
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and  then  rearrange  to  obtain 


(9) 


E[t2].  =  2 
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Ideally,  we  would  know  X  for  each  of  the  n  tissue  pieces, 
but  without  this  information  we  can  at  least  determine 
the  range  of  RDs  possible  within  this  model.  Because 
XL  >  XNI ,  assigning  the  largest  r  values  (the  lowest  blood 
flows)  in  Eq.  10  to  the  lipid  portion  maximizes  the  pre¬ 
dicted  RD,  and  assigning  the  smallest  r  values  to  the  lipid 
minimizes  it.  If  the  multiple-compartment  model  is  accu¬ 
rate,  the  distribution  of  lipid  and  blood  flows  actually 
present  in  ihe  tissue  must  lie  between  these  two  ex¬ 
tremes.  We  first  calculated  the  RDs  by  assuming  that  the 
lipid  portion  should  be  paired  with  the  higher  of  the  mea- 
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sured  blood  flows.  We  then  made  the  same  calculations 
under  the  assumption  that  lipid  should  receive  the  lower 
of  the  measured  blood  flows.  We  examined  a  third  as¬ 
sumption  that  the  lipid  is  evenly  distributed  among  the 
tissue  pieces  by  simplifying  Eq.  7  to 

El*2],-  =  2  i  +  (1  -  o))XNL]  (11) 

We  estimated  the  moments  of  the  intact  side  as  the 
weighted  averages  of  the  moments  of  the  muscle  and 
nonmuscle  tissues.  The  first  moment  for  the  intact  side  is 

E[f]23  =  (1  -  5)E [t]2  -  (5)E[t]3  (12) 

and  the  second  moment  is 

E[^]23  =  (1  -  <5)E[^]2  +  (6)E[^]3  (13) 

where  5,  the  weighting  factor  for  the  contribution  of  the 
nonmuscle  layer  to  the  detected  Xe  activity,  is  calculated 
from  the  measured  tissue  layer  thicknesses,  as  described 
in  the  APPENDIX.  Finally,  the  predicted  variance  is  equal 
to  Eft2]  -  (E[f])2,  from  which  the  RD  is  calculated  as 
( variance/mean2) 1/2. 

Statistical  Analysis 

The  number  of  exponential  terms  required  to  fit  the 
Xe  data  was  selected  according  to  reductions  in  the  resid¬ 


ual  error  by  use  of  an  F  test.  The  observed  and  predicted 
means  and  RDs  were  compared  by  calculating  95%  confi¬ 
dence  intervals. 

RESULTS 

Observed  Xe  Kinetics 

A  typical  profile  of  inspired  Xe  radioactivity  and  the 
radioactivity  detected  in  the  tissue  is  shown  in  Fig.  1.  We 
fitted  the  curves  using  an  expression  with  two  exponen¬ 
tial  terms.  This  expression  fits  the  data  well  (Fig.  1).  The 
estimates  of  the  mean  residence  times  had  coefficients  of 
variation  of  ^3.5%,  similar  to  previous  results  (12).  At¬ 
tempts  at  fitting  with  one  exponential  term  produced  visi¬ 
bly  poorer  fits,  which  were  nearly  always  rejected  as  sta¬ 
tistically  less  satisfactory  than  the  two-exponential  fits 
(P  <  0.05).  Attempts  to  fit  three  exponentials  to  the  data 
produced  statistically  significant  improvements  over  the 
two-exponential  fits  in  only  2  of  16  cases.  The  mean  resi¬ 
dence  times  estimated  from  the  three-exponential  fits 
were  similar  to  those  obtained  using  two-exponential 
terms. 

The  observed  mean  residence  times  estimated  for  the 
intact  side  were  longer  than  those  estimated  for  the  iso¬ 
lated  muscles:  35.4  vs.  23.0  min  (Table  1).  The  observed 
RDs  were  similar  for  the  two  sides:  1.61  for  intact  and 
1.76  for  isolated. 


FIG.  I.  Profiles  of  Xe  radioactivity  over  time  in  inspired  air  (Xe  input)  and  tissue  (tissue  Xe  data)  for  isolated 
muscle  in  rxpl  5  and  2-exponential  fitted  curve  (2-exp  fit).  Predicted  mean  residence  time  was  5.3  min,  observed  mean 
residence  time  was  17.2  min.  Also  shown  is  tissue  Xe  curve  predicted  by  multiple-compartment  model,  with  lower 
perfusion  rates  assigned  to  lipid  (model-predicted).  Shapes  of  curves  reflect  differences  between  observed  and  pre¬ 
dicted  relative  dispersions. 
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TABLE  1.  Observed  and  predicted  Xe  mean 
residence  times 


Expt 

No. 

Isolated 

Intact 

Observed 

Predicted 

Observed 

Predicted 

/ 

21.8±0.3 

8.8 

41.2±0.6 

30.6 

2 

39.4±2.2 

23.3 

52.0±1.7 

47.2 

3 

12.0±0.4 

9.1 

25.9±1.2 

36.8 

4 

33.0±1.3 

14.6 

60.2+8.4 

30.4 

5 

17.2±0.3 

5.3 

20.7±0.4 

7.7 

6 

21.0±0.5 

7.7 

23.1  ±0.4 

8.6 

7 

15.1+0.7 

8.8 

16.6±0.5 

13.8 

8 

24.6±0.4 

10.0 

43.6±1.1 

25.8 

Avg 

23.0 

11.0 

35.4 

25.1 

Values  are  means  ±  SD  in  minutes. 


Predicted  Xe  Kinetics 

The  average  vol/vol  fraction  of  lipid  in  the  isolated 
muscle  was  6.0%  (range  3.2-13%),  similar  to  previously 
reported  values  (4,  8),  and  the  average  in  the  tissues  of 
the  intact  side  was  13%  (range  3.4-35%).  The  corre¬ 
sponding  average  tissue-blood  partition  coefficients  esti¬ 
mated  using  Eq.  6  to  account  for  the  measured  lipid  were 
0.98  (range  0.68-1.60)  for  the  isolated  muscle  and  1.61 
(range  0.73-3.60)  for  the  intact  side.  The  average  blood 
flows  were  9.5  ml  •  min-1  •  100  g_1  for  the  isolated  muscle 
and  7.3  ml  •  min  1  •  100  g“‘  for  the  intact  side.  The  muscle 
tissue  of  the  intact  side  was  on  average  2.8  cm  thick,  and 
the  nonmuscle  tissue  was  on  average  0.4  cm  thick.  As  a 
straight  percentage,  the  nonmuscle  tissue  layers  repre¬ 
sented  just  12.5%  of  the  total  thickness  of  the  intact  side. 
However,  considering  the  locations  of  the  tissues  with 
respect  to  the  detector  and  the  geometry  of  the  field  of 
view  of  the  detector,  we  estimate  that  the  nonmuscle  tis¬ 
sues  contributed  an  average  of  25%  of  the  Xe  signal  re¬ 
corded  from  the  intact  side  (see  appendix). 

The  predicted  mean  residence  times  for  the  intact  side 
were  longer  than  those  of  the  isolated  muscles:  25.1  and 
11.0  min,  respectively  (Table  1).  With  incorporation  of 
the  perfusion  heterogeneity  alone,  the  model  predicted 
an  average  RD  of  only  1.16.  Adding  the  effect  of  lipid 
increased  the  predicted  RDs  to  2.5-5.  Pairing  the  lipid 
with  the  highest  of  the  measured  flows  produced  RDs  of 
~2.5,  whereas  pairing  the  lipid  with  the  lowest  flows  in¬ 
creased  the  predicted  RDs  to  ~5.  Distributing  the  lipid 
evenly  throughout  the  tissue  produced  intermediate  val¬ 
ues  averaging  ~3. 

Observed  vs.  Predicted  Kinetics 

The  observed  mean  residence  times  were  on  average 
twice  as  long  as  those  predicted  by  the  model  for  the 
isolated  and  the  intact  sides  (95%  confidence  interval 
1. 6-2.3).  This  can  be  appreciated  visually  in  the  plots  of 
the  mean  residence  times  shown  in  Fig.  2.  The  predicted 
mean  residence  times  fall  short  of  the  line  of  equality  in 
all  but  one  case  (Expt  3,  intact,  Table  1).  The  observed 
and  predicted  mean  residence  times  are  correlated  (r  = 
0.87  for  isolated,  r  =  0.69  for  intact;  P  <  0.05). 

Figure  3  summarizes  the  observed  RDs  and  those  pre¬ 
dicted  by  the  model.  The  observed  RD  of  1.68  is  similar  to 
previously  reported  values  (12,  16).  Incorporating  the 


measured  perfusion  heterogeneity  and  a  single  constant 
X  for  the  tissue,  ignoring  the  compositional  heterogene¬ 
ity,  the  model  predicted  an  RD  of  only  1.16,  which  is  not 
much  larger  than  the  RD  of  1  for  a  well-mixed  compart¬ 
ment.  However,  incorporating  X  values  to  reflect  the  lipid 
content  of  the  tissue  increased  the  predicted  RDs  consid¬ 
erably.  Distributing  the  lipid  evenly  among  the  heteroge¬ 
neous  perfusion  rates  produced  an  average  RD  of  3.34 
(Eq.  11).  We  estimated  the  range  of  RDs  possible  under 
this  model  by  pairing  the  lipid  with  the  higher  and  with 
the  lower  of  the  measured  blood  flows  (Eq.  10),  which 
produced  average  RDs  of  2.53  and  4.88,  respectively. 

DISCUSSION 
Mean  Residence  Time 

Lipid  alone  cannot  account  for  the  observed  residence 
times  of  Xe  in  tissue.  The  mean  residence  times  pre¬ 
dicted  by  the  present  model  were  as  long  as  physiologi- 
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fig.  2.  Plot  of  observed  (.t  axis)  vs.  predicted  (v-axis)  mean  resi¬ 
dence  times  for  isolated  (A)  and  intact  (fi)  muscles.  In  both  cases,  fitted 
regression  (dashed)  lines  fall  below  (solid)  line  of  equality,  reflecting 
longer  observed  than  predicted  Xe  residence  times.  Average  ratios  of 
observed  to  predicted  mean  residence  times  are  2.2  for  isolated  side  and 
1.7  for  intact  side. 


2132 


LIPID  AND  Xe  EXCHANGE  IN  MUSCLE 


FIG.  3.  Observed  and  predicted  Xe  relative  dispersions  (RDs),  with 
95%  confidence  intervals  shown  in  parentheses.  Average  of  observed 
RDs  (Observed)  was  1.68  (1.55-1.82).  With  incorporation  of  tissue 
blood  flow  heterogeneity  alone  (flow  only),  model  predicted  an  RD  of 
1.16  (1.11-1.21).  Including  lipid  evenly  throughout  tissue  (uniform 
lipid)  increased  predicted  RD  to  3.34  (3.10-3.58).  Pairing  lipid  with 
highest  of  measured  blood  flows  (high  flow)  predicted  an  RD  of  2.53 
(2.32-2.74),  whereas  pairing  with  lowest  flows  (low  flow)  predicted  an 
RD  of  4.88  (4.35-5.40).  As  reference,  dashed  line  shows  RD  of  1  for  a 
single  mixed  compartment. 

cally  plausible  on  the  basis  of  our  measured  tissue  lipid 
contents.  However,  our  observed  mean  residence  times 
were  twice  as  long  as  those  predicted  by  the  model.  Re¬ 
calculating  the  residence  times  with  the  A  values  held 
constant  at  published  values  of  0.68  (17)  for  muscle  and 
0.87  for  subcutaneous  tissue  (10)  gives  observed  means 
that  are  on  average  2.5  times  longer  than  the  predicted 
means. 

How  large  would  the  tissue-blood  partition  coefficients 
have  to  be  to  account  for  the  observed  mean  residence 
times  by  use  of  the  present  model?  Setting  Eq.  5  equal  to 
the  observed  mean  residence  time  and  solving  for  an 
average  A  for  the  tissue,  we  estimate  that  the  partition 
coefficients  would  have  to  be  2.2  for  the  isolated  muscle 
and  3.2  for  the  intact  side.  Both  of  these  values  are  con¬ 
siderably  larger  than  any  of  the  reported  coefficients  for 
these  tissues  (10,  17).  Alternatively,  using  Eq.  5  with  the 
average  published  A  values,  we  can  estimate  the  amounts 
of  lipid  that  would  be  required  to  explain  the  observed 
mean  residence  times:  16.5%  for  the  isolated  muscle  and 
27.5%  for  the  intact  side.  Again,  these  greatly  exceed  our 
measured  lipid  fractions  as  well  as  those  reported  in  the 
literature  (4,  8). 

Dispersion  of  Residence  Times 

The  results  shown  in  Fig.  3  suggest  that  lipid  contrib¬ 
utes  significantly  to  the  RD  of  Xe  residence  times.  The 
observed  RDs  are  larger  than  can  be  accounted  for  by 
perfusion  heterogeneity  but  are  smaller  than  the  pre¬ 


dicted  RDs  when  the  compartments  include  lipid.  What 
would  be  the  effect  of  diffusion  of  gas  between  regions? 
Computer  simulations  show  that  increasing  the  diffusion 
between  heterogeneous  features  such  as  lipid  and  nonli¬ 
pid  tissues  diminishes  the  RD  (5).  This  suggests  that  a 
strict  compartmental  representation  of  inert  gas  ex¬ 
change  is  unrealistic,  because  it  fails  to  account  for  the 
effects  of  diffusion. 

Considerations 

Because  our  results  are  not  explained  by  a  model  using 
physiologically  plausible  parameter  values,  we  must  con¬ 
sider  potential  errors  in  the  data  measurement  and  in  the 
model  itself. 

Error  in  measured  blood  flow  rate.  During  5-h  experi¬ 
ments,  do  four  blood  flow  measurements  adequately  rep¬ 
resent  the  tissue  blood  flow,  compared  with  the  hundreds 
of  data  points  of  the  Xe  curves?  If  the  tissue  blood  flow 
rates  during  the  intervals  between  microsphere  injec¬ 
tions  were  actually  half  as  large  as  those  found  by  our 
measurements,  our  biased  high  flow  averages  would  pre¬ 
dict  inappropriately  short  mean  residence  times.  Analyz¬ 
ing  the  variability  of  the  measured  flows,  we  find  that  the 
95%  confidence  intervals  are  within  one-third  of  the 
mean  flows  of  the  four  sampled  values. 

Could  smaller  pieces  of  tissue  contain  sufficiently  large 
flow  heterogeneities  to  explain  the  observed  Xe  RDs? 
Previous  measurements  (12)  on  280-mg  tissue  pieces, 
one-fourth  the  size  of  those  measured  here,  found  even 
less  flow  heterogeneity  (RD  =  1.10)  than  in  the  ~l-g 
pieces  measured  here.  Diffusional  mixing  largely 
obscures  flow  heterogeneity  influences  in  this  size 
range  (6). 

Error  in  detector -tissue  representation.  We  did  not  de¬ 
termine  the  relative  position  of  each  tissue  piece  in  the 
field  of  view  of  the  Xe  detector.  What  would  the  pre¬ 
dicted  mean  residence  time  be  if  all  the  slowly  exchang¬ 
ing  pieces  were  located  near  the  detector  and  therefore 
contributed  disproportionately  to  the  recorded  Xe  sig¬ 
nal?  Using  the  appendix  to  recalculate  the  predicted 
mean  residence  times  under  these  circumstances,  we  find 
that  the  ratio  of  observed  to  predicted  mean  residence 
times  decreases  from  2.2  to  1.8  for  the  isolated  muscle 
and  from  1.7  to  1.4  for  the  intact  side.  So  even  our  most 
conservative  consideration  of  the  lipid  and  the  detector- 
tissue  geometry  still  leaves  a  discrepancy  between  the 
observed  Xe  mean  residence  times  and  the  predictions  of 
the  model. 

Conclusions 

The  results  presented  here  demonstrate  that  a  multi¬ 
ple-compartment  model  does  not  accurately  describe  Xe 
gas  exchange  in  tissue.  Plausible  values  for  the  tissue- 
blood  partition  coefficients  and  the  lipid  contents  of  the 
tissues  do  not  explain  the  observed  means  and  RDs  of  Xe 
residence  times.  Recirculation  of  Xe  in  the  tissue  by 
countercurrent  vessel  exchange  or  other  diffusion-based 
processes  may  be  important.  Experiments  considering 
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FIG.  4.  A:  schematic  showing  cross 
section  of  detector-collimator  positioned 
above  tissue  slab.  Total  radioactivity  de¬ 
tected  from  a  point  at  ( x,q )  is  calculated 
by  integrating  detector  surface  over  (r,0) 
for  all  vectors  i>.  B:  detector  view  is  es¬ 
tablished  by  constraints  of  collimator 
edges  in  plane  of  collimator.  See  text  for 
definitions  of  abbreviations. 


gas  diffusivity  and  tissue  architecture  are  needed  to  eval¬ 
uate  these  notions. 


APPENDIX 

Model  of.  Tissue-Delector  Geometry 

To  estimate  the  fraction  of  radioactivity  contributed  by  a 
tissue  layer  of  variable  thickness  and  distance  from  the  detec¬ 
tor,  we  approximate  the  tissue  as  a  slab  with  a  uniform  density 
of  radioactivity.  A  detector  with  a  circular  surface  of  radius  R 
and  area  A  is  positioned  above  and  normal  to  the  slab  surface 
(Fig.  4).  The  detector  shaft  is  housed  ip  a  collimator  that  ex¬ 
tends  lepgth  (L)  beyond  the  detector  tip,  with  a  small  gap  (G) 
between  the  detector  and  the  collimator  interior.  The  collima¬ 
tor  tip  rests- against  the  slab  surface  at  x  =  0,  where  x  is  the 
distance  into  the  slab.  The  collimator  establishes  a  conical  view 
of  the  slib  and  allows  no  radioactivity  to  penetrate  its  walls. 
The  distance  q  extends  radially  in  the  plane  of  the  slab  normal 
to  the  cone  axis. 

A  point  at  (x,q)  establishes  a  vector  (v)  of  radioactivity  to  a 
small  area  (dA)  of  the  detector  surface  at  (r^,<t>).  The  efficiency 
of  detection  of  a  given  point  depends  on  three  factors:  the  ra¬ 
dial  location  in  the  slab  plane,  the  distance  from  the  detector, 
and  the  attenuation  of  the  radioactivity  by  the  tissue.  First,  the 
radial  location,  a  point  at  (x,q),  is  detected  with  an  efficiency 
proportional  to  cos  0,  where  cos  0  =  (x  +  L)h.  That  is,  points 
directly  in  front  of  the  detector  surface  (toward  cos  0  =  1)  are 
detected  more  efficiently  than  those  to  either  side  (toward  cos 
*72  =  0).  Second,  the  distance,  the  detected  radioactivity,  falls 
off  as  the  inverse  of  the  square  of  v,  as  if  the  point  were  on  the 
surface  of  a  sphere  of  area  4jti'2.  The  detected  activity  (dE)  as  a 
function  of  radiai  location  and  distance  is  expressed  as 


cos  0AA 

dE  =  ~^r 


(Al) 


Third,  the  tissue  attenuates  the  radiation  traveling  through  it 
(14).  This  is  expressed  as 


dE(x/J)  =  exp  (-ae’)dE  (A2) 


where  a  is  the  linear  attenuation  coefficient  (in  units  of  recipro¬ 
cal  distance)  and  v'  =  v  -  (LI cos  8)  is  the  distance  traveled 
through  the  tissue. 

The  sum  E  of  the  detected  radioactivity  at  the  detector  sur¬ 
face  from  radiation  at  any  point  (x,q)  is  then  expressed  as 


Wf 

•'o^o 


Rd  exp  (-  av')  cos  0rd 


Attv2 


drdd<*> 


(A3) 


where 


v2  =  r%  +  q2  -  2 r^q  cos  <t>  +  (x  +  L)2  ( A4 ) 


and 


rc  <  R  +  G 


(A5) 


Equation  A5  means  that  only  vectors  within  the  collimator 
view  are  counted  at  the  detector  surface. 

Equation  A3  solved  by  numerical  integration  allows  us  to 
calculate  the  radioactivity  detected  from  any  point  in  the  slab. 
Integrating  over  depth  x,  radius  q,  and  angle  6  in  the  plane  of 
the  slab  shows  that  all  planes  contribute  equally,  because  the 
enlargement  of  the  letector  view  perfectly  offsets  the  decline  in 
radioactivity  with  distance.  The  attenuation  due  to  the  tissue 
itself,  however,  causes  the  nearer  layers  to  contribute  slightly 
more  than  the  farther  ones.  All  that  is  required  to  calculate  the 
contribution  of  a  given  layer  is  information  about  its  thickness 
and  its  distance  from  the  detector  surface,  which  is  provided  by 
the  tissue  layer  measurements  described  in  METHODS.  The 
small  1-g  pieces  making  up  any  one  layer  of  muscle  or  nonmus¬ 
cle  tissue  are  considered  without  regard  to  their  relative  posi¬ 
tions  within  the  layer. 

The  fraction  of  detected  Xe  counts  (5)  coming  from  the  non¬ 
muscle  layer  of  the  intact  side  is  calculated  as  the  sum  of  the 
radioactivity  corresponding  to  the  measured  thickness  divided 
by  the  total  radioactivity  detected. 
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